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Abstract. - We propose an indirect way of studying the fluctuation-dissipation relation in 
spin-glasses that only uses available susceptibility data. It is based on a dynamic extension 
of the Parisi- Toulouse approximation and a Curie- Weiss treatment of the average magnetic 
couplings. We present the results of the analysis of several sets of experimental data obtained 
from various samples. 



Introduction. — The fluctuation-dissipation theorem (FDT) relates the response of a mag- 
netic system to the magnetization correlation function at equilibrium. In its integrated form 
FDT states that 

X(t,t w ) = ^(q d -C(t,t w )) , (1) 

where the response to an applied field h held constant from a waiting-time t w up to t and the 
correlation are defined as 

x(t,t w ) = 6{m(t))/8h\h=o , C(t,t w ) = (m(t)m(t w )} , 
and q^ is the long-time limit of C(t, t). 

(*)On leave from Service de Physique Statistique, Magnetisme et Supraconductivite, DRFMC, 
F-38054 Grenoble, France. 
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Glassy systems are out of equilibrium and FDT does not apply, as shown in particular 
in solvable models with infinite range interactions. Otherwise stated, there are effective 
temperatures (different from the bath-temperature) at play in aging systems Q. A simple 
modification of Eq. ([!]) consists in proposing that, for large t w , the susceptibility still depends 
on t and t w only through C, i.e. 

X (t,t w ) = x(C(t,t w )) (2) 

where x(C) is a system-dependent function The latter may be obtained from a plot of 
x(t,tw) against C(t,t w ) using t > t w as a parameter. This curve is known analytically for 
several mean-field spin-glasses || such as the Sherrington-Kirkpatrick (SK) model. Models 
with finite range interactions have been studied numerically by several groups |^|, ^| who 
obtained results that are qualitatively similar to the mean-field ones. These numerical results 
must however be taken with caution since the times that can be reached in simulations are 
relatively short. 

In principle, x an d C should be determined experimentally by measuring independently 
noise correlations and susceptibilities. This procedure has been recently used to investigate 
FDT violation in structural glasses ||. For spin-glasses, noise measurements along the lines 
of the early work of Ref. j?j are under way: their outcome will provide a most stringent test 
for spin-glass theories. 

In this paper, we shall assume that violations of FDT do occur below T g , and propose 
an indirect method for the determination of the x vs C curve using available experimental 
data. This construction gives a first glimpse of the form of this curve. It is based on some 
assumptions, notably a dynamic extension of the Parisi- Toulouse (PaT) approximation 
that hold for some spin-glasses. 

We start by noticing that, in all solvable spin-glass models, below T g , is a piecewise 

function ||: 



f ?{q.d-C) if q<C<q d (FDT regime), 

X[ ' 1 Xag(C) + i (q d - q) if q <C<q (aging regime), {) 

where the dynamical Edwards-Anderson parameter, q, and the minimal correlation in an 
applied magnetic field H, qo, are defined as 

q = lim lim C(t,t w ) , qo = lim C(t,t w ) , (4) 

t—t w — >OOt w — >OG t — >00 

the latter vanishing in zero applied field. In the SK model, Xag(C) is decreasing and has a 
downwards curvature. Numerical results indicate that, at least within the simulation times, 
the shape of the curve for the 3D Edwards- Anderson (EA) model is similar. 

The dynamical version of the PaT hypothesis || || (see also Ref. M) consists of the following 
two assertions: (i) is independent of T and H in the aging regime and (ii) q and qo only 

depend on T and H , respectively. We shall moreover assume that this approximation is good 
even at finite times (see Discussion). 

The near temperature-independence of x(C) in the aging regime has been checked numer- 
ically for the 4D EA model in Ref. ||. No checks are available for the 3D case. The validity 
of this approximation for the experimental systems will be discussed below. It will be seen 
that the PaT approximation allows us to estimate the C-dependence of the susceptibility using 
exclusively response results, thus circumventing the difficulties inherent to noise measurements. 

Our strategy is to use data taken under T and H conditions such that the system is at the 
limit of validity of FDT, i.e. C(t,t w ) = q. The point {q,x(l)} 1S the intersection between 
the straight part (FDT regime) and the curved part (aging regime) of x(C) (cf. Eq. ||). The 
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locus of the points obtained varying T and H spans a master curve x(C) which, by the PaT 
hypothesis, is field and temperature independent. The method of construction is explained 
below and illustrated in Fig l|. 

The susceptibility at the limit of the FDT regime corresponds to: 

X(q) = lira lim X (t,t w ) = -{q<i-q) ■ (5) 

t—t w — >oct w — >oc ± 

We have approximated this limit by using susceptibility data of three different types taken 
from the literature. 

Frequency dependent measurements. 

In ac-susceptibility measurements, a small ac-field of fixed frequency u> is applied. The 
in-phase susceptibility x'i^^w) is recorded as a function of temperature. For frequencies in 
the range oj ~ 1Hz, the long waiting-time limit tot w ^ 1 is approached within the measurement 
time. Then we can estimate the limit of zero frequency by extrapolation 



X '(0,oo) = lim km t w ) = -{q d - q(T)) . (6) 

The master curve x(C) is obtained by joining the points {C = qd — T X '(0, oo) ; X — x'(0, oo)} 
using T as parameter. 

Field cooled measurements . M FC is measured by cooling the sample in a constant magnetic 
field. Below T g , Xfc = dM FC /dH rapidly reaches an asymptotic value (see however pp|). 

In some spin- glasses like CuMn Jl2[ Xfc is nearly temperature-independent below T g , 
Xfc(T,H) ~ Xfc(H) as required by PaT. However, this does not hold in most systems for 
small fields and near T g where a cusp in Xfc appears. The x FC data may be used where both 
FDT and PaT hold: on the critical line (assuming there is at least a transient one, an issue 
discussed below). The parameter is now H and the master curve is spanned by the points 
{C = q d -T g (H) X Fc(H); x = Xfc(H)}. 

Zero-field cooled measurements . 

In the absence of a sufficient amount of published ac or FC data, we have resorted to include 
in the analysis data obtained in a zero-field cooled procedure. The sample is quenched in the 
absence of a field from above T g down to some low temperature. After a time t w (necessary for 
stabilization of the temperature) , a weak magnetic field H is applied and A/ ZFC is immediately 
measured. Under typical experimental conditions, the measurement time t is significantly 
shorter than t w . We shall thus consider that Xzfc = M ZFC /H is a good approximation to 
Eq. (§. 

Ideally, the same procedure should be repeated for each measurement temperature. In 
practice, however, the magnetization M ZFC is measured by increasing the temperature in steps 
from its initial value. Although the two methods are not strictly equivalent, we believe that the 
possible differences are of little consequence for our conclusions. Therefore we shall consider 
that the whole experimental M ZFC (T) curve yields an acceptable approximation to Eq. (^|). 
Support for this point of view is given by a comparison of both ac and zfc procedures on one 
sample (see below). 

The field-independence of x ZFC implied by PaT is in general well verified experimentally 
except for the largest fields (see, for example, Fig. 1 of the first of Refs. [0). The master 
curve is determined as in the accase. 

Before turning to the analysis of the data, we notice that Eq. ([|) cannot be applied to 
experimental data as it stands. Indeed, this equation is only valid for systems in which 
the exchange coupling averages to zero. In real spin-glasses, however, this average is, in 
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Fig. 1. - Sketch of the \ vs C plot. The thick curve represents the master curve x(C) that, within 
the PaT approximation, is temperature and field independent. The thin straight line has slope — 1/T 
(T < T g ) and represents the first line in Eq. (^|). The dashed straight line has a slope —\/T g {H) and 
joins {q d ,0) to (xFc(H),q (H)). 

Fig. 2. - x(C) plot for CuMn at 1% and 2%. Data are taken from reference 113] . The vertical axis 
is normalized by the susceptibility at the critical temperature in zero field (xo> The horizontal axis 
is normalized by q& = C. The crosses are numerical results for the 3D EA model H- The inset shows 
the inverse FC and ZFC susceptibilities as functions of temperature. 



general, finite. This is reflected in the behavior of the high temperature susceptibility that 
often obeys a Curie- Weiss law. The Curie- Weiss temperature 9 may be a sizeable fraction of 
T g . Within a mean-field approximation of the average coupling, however, Eq. (B) still holds 
for the response to the total field (applied plus internal). This amounts to the replacement 
X ~ * Xmeas/ (1 + 0/C Xmeas)) with Xmeas the measured susceptibility and C the Curie constant. 
In some metallic spin-glasses, like CuMn at low concentration of the magnetic impurity 113] , as 



well as AuFe ]13j and AgMn 1 14 , the Curie- Weiss law holds down to the transition temperature. 
This is not the case in other spin-glasses where, due to progressive clustering, the paramagnetic 
behaviour deviates from a simple Curie- Weiss law close to the transition. When the Curie- 
Weiss law holds, a plot of the inverse susceptibility as function of temperature yields the values 
of C and 9 and it is easy to show that C = q e i- 

The analysis. — We have analysed data obtained with the methods described above for 
three metallic systems: CuMn [jl2| and AuFe fi~3| f or several concentrations, AgMn at 2.6% 
[ [l4"| , and two insulating samples CdCr1.7Ino.3S4 |l5|, [l6| and Fe . 5 Mn .sTiO3 pl| . 

For CdCr1.7Ino.3S4 we used the x'(o;, t w ) data ]l6| ], and extrapolated to zero frequency 
assuming a power-law decay x'i^i 00 ) ~ x'(0, 00) + c\uj a . In the cases of CdCri. 7I.n0. 3S4 and 
Feo.5 Mno.sTiOs, we also used the FC data to check the consistency of our determination. 
The susceptibilities Xzfc were obtained for CuMn, AgMn, AuFe, CdCri. 7I.n0. 3S4 and Feo.5 
Mn . 5 TiO 3 using Xzfc ~ M ZFC /H. 
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In order to analyse the field-cooled data, it is necessary to identify the transition temperature 
in the presence of an applied field. This is done by studying the onset of irreversibilities in the 
magnetization curves: we defined T g as the temperature where M mR = M FC — M ZF c = 0. At 
high fields, when the relationship between Xfc and M FC is non-linear, we have made polynomial 
fits of M FC {H) to compute Xfc(H) = dM FC /dH. For high enough fields, M FC is independent 
of temperature supporting the PaT approximation. However, for some samples, a cusp in 
the low-field susceptibility may appear near the transition temperature and the determination 
of a T-independent Xfc becomes ambiguous. In these cases we have used the value Xfc at 
the critical temperature in order to ensure consistency with the alternative determination of 
x(C) based on M ZFC measurements. On the whole, the construction using FC data is more 
ambiguous than that using zfc or ac data. 

Our analysis is most reliable for CuMn, a system in which the Curie- Weiss law as well as 
the PaT approximation are very well verified. Figure || shows the x{C) curve determined using 
the ZFC data of Ref. (l2) for two concentrations, 1.08 % and 2.02%. There are no experimental 
points for C/qd > 0.8 that correspond to rather low temperatures. We know however that 
x(C) tends to zero as C — > qd since Xzfc{T = 0) = 0. In addition, the slope dx/dC should be 
infinite at C — qd so that q — qd only at T = 0. The validity of the hypotheses can be judged 
by the inset of Fig. || where we show the temperature dependence of the inverse susceptibiliy 
for the 1.08% compound. A Curie- Weiss law with « holds accurately for all T > T g . 
The T-independence of Xfo required by the PaT approximation is also well verified below the 
transition. The same is true for the 2.02% sample. 
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Fig. 3. - x(C) curves for CdCr1.7Ino.3S4 determined from FC, ZFC and ac susceptibilities. We show 
for comparison the ZFC data for CuMn 1.08%. 

Fig. 4. - x(C) curves for some of the samples studied. ZFC data for CuMn at 1.08% and 2.02 %; 
AuFe at 1.0%, 2.0%, 4% and 8.0%; AgMn at 2.6% and Fe .s Mn . 5 TiO 3 ; FC data for Fe . 5 Mn .sTiO 3 . 



For comparison, we also show in Fig. Q the curve x{C) for the 3D EA model, at T = 0.7(< 
T g ) and H — 0, obtained numerically in Ref. 0. The agreement between the numerical results 



6 



EUROPHYSICS LETTERS 



and the experimental data for the 1.08% sample is remarkable. It may be fortuitious, however, 
since the results for the 2.02% sample deviate from it. In fact, one must note that x(C) is not 
a universal function. For example, it depends on the details of the Hamiltonian (Heisenberg, 
Ising and, in general, the level of anisotropy) even at the mean-field level. Thus, there is no 
reason to expect universality in real systems. 

The data for CdCri. 7^0.384 |l5|, |l6) obtained using the three different techniques described 
above are shown in Fig. ^. It can be noticed that the FC and ac results are very similar. The 
ZFC data are somewhat higher (probably due to the fact that the large t w condition is not as 
well full filled as in the other cases), but the agreement with the other determinations remains 
acceptable. Notice that the x(C) curve obtained for this compound is quite different from 
that corresponding to CuMn. 

In Fig. U we collect results for all the other samples. As expected, the curves do not fall on 
a universal curve, but their shape is similar. 

Discussion. — Finally, let us clarify an important point. In several situations, such as 2D 
Ising and 3D Heisenberg |T^1 and, perhaps, 3D Ising spin glasses under a magnetic field, no true 
spin-glass phase is expected. However, for still relatively long times the system remains below 
a slowly time-dependent pseudo de Almeida-Thouless (AT) line [p8| : it ages and behaves as a 
true (out of equilibrium) glass with a non-trivial x(C) that would eventually become a straight 
line with slope —1/T. In this paper we explored the consequences of the stronger assumption 
that PaT is a good approximation below the pseudo AT-line if all quantities involved belong 
to the same epochs. 

Another important issue is the asymptotic (t w — > 00) form of the Xag(C) curve. Even if 
the system never equilibrates, the Xag(C) curve may still be a very slowly varying function of 
t w , eventually reaching a form different from that observed experimentally. We are not in a 
position to discard this possibility. 

It has been recently shown that, under certain hypotheses, the slope of the dynamic x(C)> 
for an infinite system in the large-^u, limit coincides with the static x(q) as defined by the 
probability of overlaps of configurations taken with the Gibbs measure jl9j . Since we do not 
address here the issue as to whether the AT-line and the PaT approximation survive beyond 
experimental times we cannot make any statements concerning the relation of this dynamical 
function to the corresponding equilibrium Parisi function. 

In conclusion we have presented an approximate determination of the x(C) curve char- 
acterising the deviations from equilibrium of spin-glasses through the violations of FDT, 
assuming they exist. This construction does not replace a true determination via simultaneous 
measurements of susceptibility and noise correlation, as the ones in Ref. §, 0, but it yields 
some insight into how this curve might look in reality. 

We wish to thank D. Sherrington for early discussions on this subject as well as H. Aruga- 
Katori, N. Bontemps, A. Ito, H. Takayama, A. Tobo and H. Yoshino. We also thank J. 
Hammann and M. Ocio for a critical reading of the manuscript. LFC, DRG and EV thank 
the Monsbusho contract for partial financial support. JK was partially supported by the 
'Programme Thematique Materiaux', Region Rhone- Alpes. 
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